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Cylindrical  coordinates.  Fig.  1 

Longitudinal,  tangential  and  radial  displacements  of  the  cylindrical  shell. 
Fig.  1.  Note  that  a positive  displacement  w is  inward. 

Coefficients  defining  the  shape  of  modes  of  the  shell  in  vacuo,  Eq.  (1) 

Subscript 

Radius  of  cylinder 

Velocity  of  sound  in  medium 

Velocity  of  sound  in  shell  material 

Subscript,  refers  to  damping 

r -dependent  function  in  expression  for  potential,  Eq.  (8) 

Distribution  of  force  P in  z -direction 
Fourier  Inverse  of  g(z) 

Hankel  Function  of  second  kind 
Thickness  of  shell 

Bessel  Function  of  first  kind 

Length  of  longitudinal  half -wave  of  mode  of  vibration 

Spacing  of  ring  stiffeners  of  shell 

Mass  of  shell  per  unit  area 

Virtual  mass  of  entrained  medium 

Coefficient,  see  Eq.  (6) 

Number  of  circumferential  waves  of  mode  of  vibration 
External  radial  force,  see  Eqs.  (3),  (11) 

Pressure  in  medium  due  to  force  P and  P,  respectively 
Subscript,  refers  to  resonance 
P'orce  exerted  by  diaphragms 
Defines  spacing  of  diaphragms.  Fig.  9 


^Additional  symbols  in  the  appendices  are  defined  as  they  occur. 
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Time 


Bessel  Functions  of  second  kind 
Coefficient,  Eq.  (7) 

Coefficient,  see  Eq.  (9) 

Coefficient,  see  Eq.  (10)  and  Table  A 
Damping  coefficient 
Dirac's  Delta  Function 
Decay  constant,  Eq.  (35) 

Variable  of  integration  related  to  L,  Eq.  (2) 
Mass  density  of  medium 
Frequency  of  vibration 
Frequency  of  cylinder  in  vacuo 


Synopsis 


The  purpose  of  the  paper  is  the  determination  of  the  sound  field  due  to  a source 
inside  a cylindrical  shell  which  is  submerged  in  an  acoustic  medium.  Using  results 
from  a previous  publication  (Ref.  1)  it  is  possible  to  express  the  radiation  by  a Fourier 
Integral;  asymptotic  integration  for  points  at  large  distances  furnishes  surprisingly 
simple  results. 

While  the  fundamental  theory  applies  for  infinitely  long  unstiffened  shells,  it  is 
extended  with  approximations  first  to  stiffened  shells,  and  as  a final  step  in  a qualita- 
tive manner  to  finite  shells. 

Particular  attention  has  been  given  to  the  question  of  resonant  frequencies  i.e. 
those  at  which  unusually  high  pressures  may  occur.  For  this  purpose  it  became  neces- 
sary to  extend  the  theory  to  Include  the  effects  of  structural  damping;  interesting  con- 
clusions concerning  the  radiation  were  obtained  by  considering  the  interaction  between 
structural  and  radiation  damping. 

Numerical  results  giving  typical  pressure  distributions  at  large  distances  were 
obtained  for  a lightly  stiffened  steel  shell  in  water.  It  is  found  that  such  a shell  has 
resonant  frequencies  with  pronounced  peaks  of  the  sound  pressure. 

A detailed  study  on  the  effect  of  rigid  diaphragms  on  the  resonant  cases  was 
made.  It  was  found  that  the  addition  of  dia{^ragms  reduces  the  resonant  effects;  for 
the  steel  shell  selected  as  example  diaphragms  reduce  resonant  effects  to  an  insignifi- 
cant level. 
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1.  Previous  Results  on  Forced  Vibrations  of  Submerged  Unstiffened  Shells 

The  response  of  an  infinitely  long,  thin  cylindrical  shell  submerged  in  an  acous- 
tic medium  due  to  sinusoidally  distributed  periodic  radial  forces  was  obtained  in  Ref.  1, 
using  the  modes  of  free  vibrations  of  the  shell  in  vacuo  as  generalized  coordinates.  For 
each  pair  of  parameters  n and  4 three  frequencies  Wj^  exist,  k = 1,  2.  3.  the  correspond- 
ing modes  having  longitudinal,  tangential  and  radial  displacements,  respectively,  Fig.  1, 

u = ^ cos  n G cos  4 z 


V = 


w 


^ sin  n 6 sin  4 z 
cos  n G sin  4 z 


(1) 


where  U||,  V||  and  W||  are  numerical  values  which  for  certain  ranges  of  n and  4 tabu- 
lated in  Ref.  2.  The  parameter  4 defines  the  half  wave  length  L used  in  Refs.  1 and  2. 


L = 


(2) 


Let  a space  and  time  periodic  external  radial  force,  act  on  the  submerged  shell, 

P(t,  G,  z)  * cos  n G e'^*  e'”*  (3) 

where  the  complex  representation  of  P used  here  differs  from,  but  is  equivalent  to  the 
one  in  Ref.  1.  The  pressure  p in  the  fluid  at  a point  having  the  coordinates  r.  G and  z is 


F(r)  cos  n G e'^'e*^^ 
P " F(i)  — ; — ZH 


(4) 


and  the  radial  displacement  w of  the  shell  at  a point  having  the  coordinates  G and  z is 

co;nGe'4«e»“^ 


w 


’ft] 


(5) 


a is  the  radius  of  the  cylinder,  and  p and  m are  the  mass  density  of  the  acoustic  me- 
dium, and  of  the  shell  per  unit  of  area,  respectively.  N is  a function  of  the  elastic 
properties  of  the  shell 


N = 


o,n* 


- n*  «{ 


n* 


(6) 


N depends  only  on  the  three  frequencies  i>>|^  of  the  shell  in  vacuo  having  half  wave  length 
L = |ir/g|  . The  coefficients  depend  on  the  shape  of  the  modes 


WJ 


W. 


r (k  • I 2.  3) 


(7) 


and  are  also  tabulated  in  Ref.  2.  If  n ■ 0 one  of  the  coefficients  vanishes  and  Eq.  (6) 
has  only  two  terms. 

The  function  F(r)  in  Eq.  (4)  arises  in  the  solution  of  the  wave  equation  for  the 
acoustic  medium,  and  is  defined  by  Eq.  (42)  of  Ref.  1 in  different  ways  depending  on  the 
value  of  n.  This  is  inconvenient  for  the  present  purpose  and  can  be  avoided  by  a modi- 
fication. Noting  that,  for  a real  x,  the  Hankel  Functions  of  the  first  and  second  kind  are 
related  H,{'*  (ix)  = (-if**  (-ix)  one  obtains  for  the  entire  range 


h|,*‘  (pa) 


where  the  sign  of  the  square  root  must  be  selected  such  that* 


[Re(p)>0.  Im(p)<0] 


(«) 


(9) 


Finally,  the  term  A occuring  in  Eq.  (4)  and  (S)  is  also  a function  of  p. 


a « A 


(P*) 


(Pa) 


(10) 


Numerical  values  of  A for  real  values  of  the  argiunent  will  be  required  later.  Such 
values  are  tabulated  in  Table  A. 


11.  Unstiffened  Shells  Acted  on  by  Forces  Close  to  the  Origin  a » 0 

The  solution  (4)  will  be  utilized  in  this  paper  to  find  the  pressure  field  due  to 
radial  periodic  forces  P(t,  0,  z)  which  are  sinusoidally  distributed  in  the  circumferential 
direction,  but  non-periodic  in  the  s -direction 


*Note  that  p as  defined  here  corresponds  to  p in  Ref.  1. 
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P(t,  9,  z)  > cos  n 0 g(z)e 


(11) 


int 


It  is  intended  to  consider  forces  acting  close  to  the  origin  z ■ 0 such  that 

I |i(*)|  d* 

-« 

is  finite.  The  restriction  on  the  integral  over  g(z)  excludes  cases  of  forces  periodically 
spaced  in  the  z -direction. 

Using  the  Fourier  inverse  of  g(z) 


0(6)  ■ 


g(z)e*^*  dz 


(12) 


the  loading  function  P(t,  0.  z)  can  be  expressed  by  the  Integral 


P(t.  0.  z)  = cos  n 9 e“” 


^ G(e)e***de 

-QO 


(13) 


This  equation  expresses  P as  an  integral  over  a spectrum  of  forces  of  the  type  of  Eq. 
(3),  and  the  fluid  pressure  p due  to  the  forces  P can  therefore  be  expressed  by  an  inte- 
gral over  the  spectnun  of  pressures,  Eq.  (4);  substituting  Eq.  (8)  the  pressure  p at  a 
point  having  the  coordinates  r,  9,  z is  given  by  the  integral 


P * 


cos  n 9 e***^ 


ao 


1 


h1*’  (pr)  G(e)e‘»*  de 

iw  [‘  - pSNljii!)] 


(14) 


It  is  noted  that  the  integrand  has  branch  points  at  | + G/c  . On  examination  it  is  found 

that  the  restriction  on  p . Eq.  (9).  requires  the  path  to  pass  above  the  branch  point  at 
G/c  and  below  at  - G/c . It  is  noted  further  that  poles  may  occur  on  the  real  axis,  i.e. 
on  the  path  of  integration.  The  required  value  of  the  integral  is  the  principal  one.  i.e.  it 
is  the  limit  for  • 0 of  the  sum  of  the  integrals  taken  over  the  sections  of  the  real  axis 

obtained  by  excluding  those  portions  of  the  real  axis  which  are  less  than  a small  distance 
f from  the  poles.  The  poles  occur  wherever  the  denominator  bracket  in  Eq.  (14)  vanishes; 
by  comparison  with  Ref.  1,  Eq.  (28),  it  is  seen  that  the  bracket  is  identical  with  the  fre- 
quency equation  for  free  vibrations,  and  poles  occur  at  points  where  the  forcing  frequency 
G coincides  with  a frequency  of  free  vibration  of  the  shell  having  a half  wave  length 
L » I v/i  I . It  has  been  shown  in  Ref.  1 that  free  vibrations  (G  and  ( real)  occur  only  if 
G4  vc/L  > I a(  I and  poles  are  therefore  only  possible  for  | i | > G/a.  The  denominator 
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i(!^ 


1 


/ 


! 


r 


bracket  being  an  even  function  poles  must  occur  in  pairs,  at  ± (j ; there  may  be  none, 
or  up  to  three  pairs  depending  on  the  value  of  n.  A typical  path  of  integration  is  shown 
in  Fig.  2. 

The  integrand  of  Eq.  (14)  is  a highly  complicated  function  of  the  variable  of  inte- 
gration 4.  and  an  exact  evaluation  of  the  formal  solution  (14)  is  out  of  the  question. 
However,  it  is  possible  to  evaluate  the  integral  for  large  values  of  r in  closed  form  as 
shown  in  Appendix  A. 

The  asymptotic  value  of  the  pressure  is 


P = exp 


2 cos  n 6 e*®*G(60) 

-paN(L“)a(a^)] 


where 


(15) 


zn 


c Vr*  + z* 

L .wcVr*f  z* 
0 (iz 


(15a) 


(15b) 


(15c) 


The  term  N(L^)  occuring  in  Eq.  (15)  is  the  value  of  N according  to  Eq.  (6)  for  the  half 
wave  length  L^. 

Due  to  its  derivation  Eq.  (15)  applies  only  for  sufficiently  large  values  of  the  ra- 
dius r.  To  get  at  least  a crude  approximation  r must  be  large  enough  to  satisfy  the  two 
conditions  (see  Appendix  A) 


- -I.*- — > (ifn<6) 


(15d) 


Ehccept  for  very  low  frequencies,  or  large  ratios  s/r,  the  first  condition  controls. 

An  alternative  asymptotic  evaluation  for  large  values  of  s is  feasible;  it  would 
essentially  give  information  on  the  conduction  of  sound  in  the  direction  of  the  axis  of  the 
shell.  This  paper  being  aimed  at  other  questions,  the  matter  is  not  carried  through. 

If  the  applied  load  acts  in  the  vicinity  of  z ■ 0 it  is  to  be  expected  for  physical  rea- 
sons that  the  pressure  at  large  distances  will  not  depend  greatly  on  the  distribution  of  the 
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applied  load  P(t,  0,  z)  in  the  z -direction  and  we  consider  therefore  only  the  effect  of  a 
line  load  acting  on  the  circle,  z ■ 0, 

P(t,  0.  z)  « 4(z)cos  n 0 e‘“^  (16) 

where  6(z)  is  Dirac's  6-function:  Eq.  (12)  gives 

G(6)  - ^ (17) 


The  peak  values  of  the  pressure.  | p |.  for  the  case  of  a line  load  according  to 
Eq.  (16)  becomes 


-I  _ cos  n 0 

1 

/ 

ir  Vr*  + z* 

hJ,®  (aPo) 

1 _ 2m 

paN(Lp)A(af^) 

For  z > 0 Eq.  (18)  simplifies  further,  as  * 0 and  > « become  independent  of  the 
frequency  n, 


- I cos  n Q 

Ph-TF— 


H 


«> 


so 

c 


I paN(«o)a 


(19) 


111.  Application  to  Shells  with  Closely  Spaced  Ring  Stiffeners 

It  was  pointed  out  in  Ref.  1,  p.  15,  that  the  results  for  the  unstiffened  shell,  quoted 
in  Sec.  I,  are  applicable,  approximately,  to  shells  with  closely  spaced  ring  stiffeners. 
Fig.  3,  if  the  frequencies  of  the  stiffened  shell  in  vacuo  are  used  as  values  occuring 
in  the  expression  for  N,  and  if  the  mass  of  the  stiffeners  is  averaged  over  the  surface 
of  the  shell  and  included  in  the  mass  term  m.  The  approximation  can  be  expected  to  be 
good  provided  the  length  L of  the  half  wave  is  several  times  the  stiffener  spacing  i,  and 
provided  the  forcing  frequency  (I  is  sufficiently  low,  such  that  the  local  displacement  of 
the  shell  between  the  stiffeners  is  unimportant. 

The  results  of  the  present  paper  were  obtained  by  integration  of  contributions 
coming  from  all  wave  lengths  L > | ir/|  |,  including  the  range  L.  < t,  and  this  appears  at 
first  to  preclude  the  application  to  stiffened  shells.  However,  in  the  range  where  the 
integrand  is  a poor  approximation. 


P “ 
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is  imaginary  if  n is  sufficiently  small,  i.e.  if  0 < tc//.  If  p is  imaginary  the  term 


lim 

r -»oo 


Hn”  (Pr) 


j/a^^r-a)lpl 


(20) 


in  the  integrand  of  Eq.  (14)  decays  exponentially  with  increase  in  r,  and  the  range  in 
question  gives  no  contribution  to  the  pressure  p at  large  distances.  (Physically  speak- 
ing this  is  so  because  these  components  of  the  shell  motion  do  not  radiate  energy,  see 
the  discussion  of  forced  vibrations  if  f)  < vcA^  Ref.  1.  p.  12.) 

It  is  therefore,  concluded  that  Eq.  (15)  is  applicable  as  long  as  R < me! pro- 
vided L is  small  enough  such  that  the  local  displacements  of  the  shell  between  the  stif- 
feners are  unimportant. 

Further  consideration  indicates  that  Eq.  (15)  may  even  be  applicable  for  larger 
values  of  R.  but  only  for  locations  where  z/r  is  small.  The  essence  of  the  process  of 
asymptotic  integration  which  leads  to  Eq.  (15)  is  the  fact  that  the  pressure  at  any  point 
is  essentially  affected  by  the  radiation  of  a component  of  the  shell  motion  of  half  wave 
length  Lq.  Eq.  (15c).  While  no  proof  is  available  one  can  assume  that  the  present  re- 
sults are  meaningful  for  a stiffened  shell  provided  > / , i.e.  if 


R < 


wc 


(21a) 


The  additional  requirement  that  the  local  displacements  of  the  shell  between  the 
stiffeners  are  unimportant  will  be  satisfied  for  a value  n only  if  the  respective  lowest 
frequency  (L  > «)  of  the  shell  in  vacuo  for  L > «is  smaller  than  the  "local”  frequency 
of  the  shell  if  the  stiffeners  were  rigid.  Approximating  the  latter  by  the  frequency  of  a 
clamped  plate  of  span  L,  one  obtains 


he, 

w,(L -«)«6.2-^  (21b) 

where  c,  is  the  velocity  of  sound  in  the  shell  material.  It  is  further  necessary  that  the 
forcing  frequency  R is  smaller  than  the  local  frequency  of  tiie  shell  in  tiie  sulmaerged 
state. 


he,  f M yl/2 

R « 6.2  -j|-  fs  j (21c) 

The  last  term  is  based  on  the  approximate  expression  for  the  virtual  mass  m,  • p//e  . 

In  order  to  apply  Eq.  (15)  in  the  suggested  manner  to  ring  stiffened  shells  of  the 
type  shown  in  Fig.  3,  it  is  necessary  to  find  the  frequencies  of  such  shells  in  vacuo. 
Suitable  approximate  methods  for  this  purpose  will  be  presented  elsewhere. 


The  approximate  treatment  of  stiffened  shells  in  the  manner  indicated  is  tanta- 
mount to  considering  the  stiffened  shell  as  an  orthotropic  unstiffened  shell.  As  n be- 
comes larger  the  value  u,(L  > «)  increases  rapidly  and  the  limitation  Eq.  (21b)  indicates 
that  the  approach  can  only  be  used  for  low  values  of  n.  For  high  values  of  n the  rings 
are  relatively  so  rigid  that  they  may  be  considered  as  absolutely  rigid;  this  fact  is 
utilized  for  an  approximate  treatment  of  resonant  cases  for  large  n.  see  Secs.  Vie 
and  VHc. 


IV.  Discussion  of  the  Results  for  Infinitely  Long  Shells 


a.  General 

The  most  remarkable  fact  is  that  the  term  N(l^)  in  Eq.  (15)  expresses  the  entire 
effect  of  the  elastic  behavior  of  the  shell  on  the  pressure  p at  a point  A,  and  more  gen- 
erally on  the  energy  radiated  in  the  direction  OA.  Fig.  4.  The  term  N depends  in  turn 
only  on  the  three  frequencies  and  the  corresponding  modes  in  vacuo  of  the  one  half 
wave  length  = | | , the  three  coefficients  being  functions  of  the  shapes  of  the 

respective  mode.  This  half  wave  length  L^is  a function  of  the  ratio  r/z  and  of  the  fre- 
quency n,  but  not  of  the  circumferential  wave  number  n.  For  a given  frequency  n the 
radiation  in  any  direction  OA  is  only  affected  by  the  vibrations  of  the  shell  of  half  wave 
length  L^. 

It  is  also  noted  that  the  pressure  decays  as  the  inverse  of  the  distance  from  the 
origin.  Vr*  + z*  . This  decay  is  similar  to  the  one  found  in  the  case  of  a pulsating 
force  acting  on  a spherical  shell.  The  decay  differs  fundamentally  from  the  one  found 
in  Ref.  1 for  a cylindrical  shell  acted  on  by  forces  which  are  space  periodic  (as  sin  i z) 
in  the  direction  of  the  longitudinal  axis  z ; for  such  loads  the  decay  was  exponential  and 
very  fast  in  the  non -radiating  low  frequency  range,  but  as  1/)^  for  high  frequencies. 
The  difference  can  be  attributed  to  the  fact  that  in  the  latter  case  energy  is  radiated  es- 
sentially only  in  two  dimensions,  while  the  presently  treated  case  is  three  dimensional. 

Eq.  (18)  was  used  to  plot  typical  pressure  curves  for  n = 0,  1 and  2 for  line  loads, 
Eq.  (16),  acting  on  a lightly  stiffened*  thin  steel  shell  in  water.  Fig.  5 shows  the  pres- 
sure p versus  the  angle  arc  tg  r/z  . Due  to  the  fact  that  r must  be  larger  than  a certain 
limit,  Eq.  (15d),  the  curve  never  applies  for  directions  parallel  to  the  axis  of  the  shell. 
The  pressure  p and  the  radiation  vanish  for  certain  directions.  The  points  p > 0 can  be 
traced  to  the  fact  that  the  term  N(Lg)  in  Eq.  (18)  vanishes.  The  case  N = 0 is  discussed 
in  Ref.  1;  it  corresponds  to  a vibration  damper  situation  in  which  the  shell  vibrates 


*The  rigidity  of  the  stiffeners  selected  was  such  that  the  fundamental  frequency  for 
n > 2 for  the  mode  of  half  wave  length  L ■ i.e.  when  the  shell  vibrates  like  a ring 
in  vacuo,  was  0.16  c/a.  The  forcing  frequency  was  n ■ 0.8  c/a. 
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longitudinally  end  tangentially,  while  the  radial  diaplacementa  w vanish.  The  number 
of  points  for  which  N ■ 0 is  one,  or  less  for  n ■ 0,  and  two  or  less,  for  n > 1,  depend- 
ing on  the  frequency.  For  loads  which  are  not  line  loads,  where  Eq.  (17)  does  not  apply, 
additional  points  p ■ 0 not  included  in  the  above  statement  occur  wnerever  G((q)  ■ 0. 

b.  Resonance 

It  is  of  considerable  interest  to  examine  the  results  for  possible  resonance  at 
special  frequencies.  The  following  discussion  applies  to  steel  shells  in  water  and  may 
require  modification  for  other  media.  The  crucial  term  to  consider  is  the  expression 
in  the  bracket  in  the  denominator  of  Eqs.  (15)  and  (18).  One  can  see  immediately  that 
this  expression  can  never  vanish,  because  the  term  h(a8^ ) is  complex  for  all  admissible 
values  of  aPo.  while  all  other  terms  are  real.  However,  in  the  range  where  the  imagi- 
nary part  of  A is  very  small  there  is  a possibility  for  very  large  magnification,  if  the 
real  part  of  the  value  of  the  bracket  should  vanish  in  the  same  range.  Table  A indicates 
that  for  n y 0,  ImA  is  quite  small  for  low  values  of  the  argument,*  n/2.  Eq.  (15b) 

indicates  that  such  low  values  correspond  either  to  low  values  of  n.  n/2.  or  to 

large  values  of  s/r;  in  either  case  large  magnification  indicating  resonance  will  occur 
if  the  physical  properties  of  the  shell  and  medium  are  such  that  the  real  part  of  the  value 
of  the  bracket  vanishes  or  becomes  very  small.  The  approximation  Re  A ^ 2/n,  (good 
for  n/2)  leads  to  the  condition 

(nM)  (22) 

If  for  a given  low  frequency.  cn/2a,  a range  of  lengths  Lq  > vc/n  can  be  found  which 
satisfies  Eq.  (22)  there  will  be  very  large  radiation;  usually  Eq.  (22)  will  be  satisfied, 
if  at  all.  for  any  length  > «c/n.  indicating  large  radiation  in  all  directions. 

It  is  noted  that  Eq.  (22)  is  actually  the  frequency  equation  for  the  submerged  shell 
(Ref.  1,  Eq.  (28) ) if  A — 2/n  is  substituted.  The  real  roots  of  this  frequency  equation 
all  correspond  to  non-radiating  lengths  L^<  «c/Rand  cause  no  pressure  for  large  radii 
r.  If  the  frequency  equation  has  complex  roots  R with  small  imaginary  parts.  Eq.  (22), 
being  an  approximation  of  the  frequency  equation,  will  be  satisfied.  Any  large  magni- 
fication found  can.  therefore,  be  ascribed  to  the  proximity  of  a subtly  damped  natural 
frequency  of  the  shell. 

c.  Steel  Shells  in  Water 

For  the  case  of  thin  steel  shells  in  water  the  criterion  Eq.  (22)  indicates  resonant 
cases  for  n h 2 , but  not  for  n ■ 1 (unless  die  mass  of  die  shell  is  about  twice  or  more 


*The  limits  for  1mA  being  small  vary  with  nt  for  1 < n < 6 the  value  n/2  is  reasonable, 
for  larger  n die  range  extends  and  for  6 < n < 15  the  limit  is  2n/3 . 
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than  the  mass  of  the  displaced  water).  For  n > 0 no  resonance  with  high  magnification 
is  possible  as  was  pointed  out  before.  Figs.  6a.  b show  the  pressure  p for  z ■ 0 as  func- 
tion of  the  forcing  frequency  n for  the  shell  used  for  Fig.  5.  As  expected,  very  narrow 
but  extremely  high  peaks  occur  in  the  response  for  n > 2.  Due  to  the  fact  that  for  low 
frequencies  the  argximents  a nk  of  the  Hankel  Functions  and  of  the  function  A are  quite 
small  it  was  found  necessary  and  convenient,  to  use  in  this  range  approximate  values 
of  these  functions  for  small  arguments  rather  than  interpolate  the  tabulated  values. 

Another  aspect  of  the  resonant  cases  is  considered  in  the  next  section. 

Fig.  6a  demonstrates  an  important  trend  concerning  the  magnitude  of  the  pressure 
p as  function  of  the  circumferential  wave  number  n.  It  is  apparent,  that  for  very  low  fre- 
quencies the  pressure  for  n > 1 is  dominant,  excepting  resonant  peaks  for  n > 2;  it  can  be 
shown  easily  that  for  small  arguments  afl/c  the  pressure  p is  a linear  function  if  n > 1. 
but  a quadratic  one  for  n » 0.  and  of  the  n^  power  for  the  non -resonant  ranges  n > 2 . 

While  Figs.  6a.  b only  contain  cases  n < 5.  it  can  be  shown  for  all  values  of  n 
that  the  response  p for  low  values  of  n.  up  to  n « 2a  Clk  is  several  magnitudes  larger 
than  for  n ^ 2a  (I)b.  excluding  points  of  resonance.  If  a force  is  considered  which  is  not 
sinusoidally  distributed  over  the  circumference,  it  can  be  expanded  in  a Fourier  series; 
the  pressure  is  then  also  obtained  in  a Fourier  series,  in  which  only  a limited  number 
of  terms  will  be  significant.  Fig.  7 shows  the  pressure  in  the  plane  z 0 as  function  of 
the  angle  0 for  a concentrated  radial  unit  load  acting  at  0 < 0 ; only  three  terms  of  the 
Fourier  series,  n.$  2,  had  to  be  used  for  n - 0.8  a/c.  Fig.  7 shows  a very  unexpected 
uniformity  of  radiation  in  all  directions.  Except  for  frequencies  where  resonance  occurs 
for  one  of  the  components,  this  behavior  seems  representative  for  aflA:  ^0.5.  For 
smaller  frequencies  the  n - 1 component  is  dominant  and  the  radiation  must  vary  as 
cos  0,  giving  a different  picture. 

d.  Limitations  Due  to  Stiffener  Spacing 

The  results  shown  in  Figs.  5-7  apply  strictly  only  to  an  orthothropic  shell,  and 
only  approximately  to  a ring  stiffened  one,  see  Sec.  III.  To  find  the  limits  of  applica- 
bility to  stiffened  shells  assume  the  stiffener  spacing  ^ ■ 0.25a,  shell  thickness  h > aAOO, 
the  ratio  c,/c  > 3.5,  and  examine  the  limitations  Eqs.  (21a,  b.  c). 

Of  the  two  limits  for  the  frequency  (21a,  c)  the  first  one  gives  af)/c  < 12.5,  the 
second  one,  which  controls  in  this  case,  ank  « 2.7 ; this  indicates  that  the  portion  of 
Fig  6b  for  larger  frequencies  does  not  apply  for  the  selected  shell. 

The  limitation  Eq.  (21b)  is  of  a different  nature.  It  requires  the  frequency  of  the 
stiffening  rings  in  vacuo  for  the  particular  value  n to  be  below  the  limit  w,(L  < «)  «3.5  c/a. 
As  these  frequencies  increase  with  the  number  n,  this  limitation  excludes  hi|^  values  of 
n.  The  ring  stiffeners  were  selected  such  that  for  n ■ 2 the  value  ut(L  ■ ■ 0.16  c/a; 

the  bending  frequencies  of  a ring  increase  as  n(n*  - 1)/Vn*+  1;  the  frequency  for  n * 7 
being  w^(L  > «i)  2.8  c/a,  the  analysis  does  not  apply  for  n > 6. 

It  follows  ttiat  the  results  shown  in  Figs.  5 -7  are  valid  also  for  stiffened  shells. 
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except  for  the  high  frequency  part  of  Fig.  6b.  It  is  of  interest  to  state  that  it  seems 
possible  to  extend  the  theory,  and  eliminate  at  least  the  limitations  Eqs.  (21b  and  c).  It 
will  be  seen  in  Sec.  Vllb  that  even  the  present  theory  is  applicable  in  the  resonant  cases 
in  the  range  where  Eqs.  (21b,  c)  are  violated. 


V . Peak  Pressures  at  Resonance,  and  Effect  of  Structural  Damping 


The  largest  resonant  value  of  p for  a given  frequency  (1  will  occur  approximately 
in  the  location  z/r  for  which  the  value  (15a),  is  such  that  the  real  part  of  the 

denominator  bracket  in  Eq.  18  vanishes.  Let  this  value  ^ (q  ” < fl/c-  Consider  the 

denominator  D of  Eq.  (18). 


D 


(2] 

where  H„ 


and  a are  functions  of 


(23) 


Using  Eq.  (10)  and  the  relation  - il^  , where  J„  and  are  Bessel  Functions  of 

the  first  and  second  kind, 


D = J„  - 


afLm 


(24) 


In  the  resonant  range  the  functions  Y are  very  much  larger  than  the  J.  and  the  vanish 
ing  of  the  last  term  causes  resonance.  Using  this  condition. 


"JfT 


Y^«0. 


to  eliminate  N in  the  first  term  of  Eq.  (24)  gives 

^ w-vx.  2 

° — V — 


(25) 


where  the  relation  on  the  Wronskian  (Ref.  3,  p.  76)  was  used.  The  pressure  p is  finally 


cos  n ^ 


(26) 


As  the  values  Y'  are  very  large  if  the  argument  is  much  smaller  than  the  order 
n.  this  pressure  is  much  higher  than  the  one  at  the  same  point  at  a slightly  different 
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frequency ; the  order  of  magnitude  of  the  latter  can  be  obtained  by  noting  that  the  bracket 
In  Eq.  (23)  In  this  range  is  of  the  order  of  unity.  With  H**’  ~ Y. 


cos  n 6 


^ 1 


(27) 


The  ratio  between  the  values  of  Eq.  (26)  and  (27),  i.e.  the  magnification,  can  be 
enormous  if  ap,,  « n.  This  fact  makes  it  necessary  to  consider  the  effect  of  structural 
damping.  It  can  be  introduced  approximately*  in  the  previous  computations  by  replac- 
ing the  term  1^  wherever  it  occurs  by  (1  U)/H,  The  effect  of  damping  in  general  is 

found  to  be  insignificant,  except  in  the  resonant  cases.  Proceeding  from  Eq.  (24)  to  find 

the  new  value  of  p , one  obtains 
*0 


W r*  + r* 
cos  n 6 


~r~ 


(28) 


where  small  terms  Jj,  were  neglected  versus  Y^  terms.  For  not  too  large  arguments 
Y„Yp*  will  be  so  large  that  the  term  2/v  is  insignificant  and 


(29) 


There  is  still  a magnification,  1/6  , but  as  1/Y„  is  very  small  the  actual  pressure  level 
becomes  unimportant  compared  to  those  for  other  values  n.  However,  if  the  resonance 
occurs  for  a value  ap^  close  to  n,  but  still  such  that  Ima  is  small,  the  situation  is  dif- 
ferent and  high  pressures  are  possible  despite  damping. 

An  estimate  for  the  maximum  value  of  ^ is  obtainable.  For  n » ap 


is  a reasonably  good  approximation;  substituting,  Eq.  (28)  has  a maximum  for 


* The  correct  way  to  Introduce  structural  damping  is  to  introduce  damping  terms 
2A  w^qj^  in  the  equations  of  motion  (Rsf.  1.  Eq.  (34) ),  where  6||  > 0 is  the  ratio  of  critical 
damping.  This  leads  to  a new  expression  replacing  the  terna  N, 


N.  - 5^ 


a^n 


As  n,  and  w|,*are  necessarily  positive  it  is  easy  to  show  that  1/Np  must  have  a posi- 
tive real  part  and  (me  can  use  the  approximation  l/Np-**  (1  + 16)/N. 


-12- 


Due  to  the  approximation  this  is  an  upper  limit  for  the  actual  value  of  p^^  . 


irVr*  + z* 

cos  n B 


max 


(30) 


The  maximum  values  of  |p^  | increase  for  a given  6 with  Vn.  In  order  to  Judge 
their  importance  they  should  be  compared  not  only  to  the  values  of  Eq.  (27)  for  the  same 
number  n,  but  also  to  the  values  of  the  pressure  p for  the  frequency  D but  different 
values  of  n.  This  changes  the  picture  because  for  a given  frequency  aD/c  the  modes 
n^  an/c  have  much  larger  values  p than  those  where  n > aD/c,  particularly  if  aD/c  is 
large. 

The  level  of  the  pressures  for  high  frequencies  in  the  lower  modes  can  be  esti- 
mated by  the  following  approximation:  if  the  frequency  D is  sufficiently  large,  a(l/c  » n, 
the  effect  of  the  elastic  resistance  of  the  structure  vanishes  and  essentially  only  the  ef- 
fect of  the  mass  of  the  shell  remains;  this  expresses  itself  in  the  previous  equations  by 
N - 1.  Using  the  asymptotic  expression  Eq.  (A -3)  gives  the  comparitive  level  of  non- 
resonant pressures 


n-.lPl-fsrVli 

The  peak  value,  Eq.  (30),  occurs  in  the  vicinity  of  ap  > n,  and  it  is  seen  imme- 
diately that  the  pressure  given  by  Eq.  (30),  at  least  for  sufficiently  large  n,  will  exceed 
the  value  of  Eq.  (31). 

Before  discussing  the  results  further,  attention  is  drawn  to  the  fact  that  Eqs.  (26), 
(28)  and  (30)  give  the  pressure  p for  any  frequency  D for  which  a resonance  occurs. 

The  value  of  the  frequency,  if  any,  for  which  this  happens  does  not  follow  from  these 
equations,  but  must  be  found  from  Eqs.  (18,  19).  This  means,  even  if  a resonance  occurs 
for  a given  shell  it  may  be  at  a point  where  Eq.  (28)  gives  a very  much  lower  value  than 
the  limit,  Eq.  (30). 

To  illustrate  the  effect  of  structural  damping  on  the  pressures  p at  resonance. 
Figs.  8a,  b show  the  values  of  p at  resonance  without  damping,  and  with  structural 
damping  6 * 0.01,  for  n - 2 and  10.  For  comparison  the  level  of  nonresonant  pressures 
for  the  steel  shell  considered  in  Sec.  IV  is  also  shown.  In  spite  of  the  low  structural 
damping  value  selected  the  change  for  low  values  aPn  is  revolutionary.  While  the  un- 
damped value  goes  to  « if  aPn  -•  0,  the  damped  values  vanish  if  aPj,  -*■  0.  As  ap^  approaches 
n the  two  curves  merge  and  damping  becomes  unimportant.  The  damped  curve  has  a 
peak,  approximately  given  by  Eq.  (30). 

The  enormous  effect  of  low  structural  damping  for  aP||  « n can  be  interpreted 
physically.  As  aPj,  -»  0,  radiation  damping  goes  to  zero,  while  the  structural  damping 
remains  finite.  As  the  energy  infnit  will  divide  in  the  ratio  of  the  damping  terms,  all 
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energy  is  absorbed  by  the  structure  and  the  radiation  pressure  p must  vanish.* 

The  large  effect  of  structural  damping  for  small  values  of  aPn  eliminates  the 
effect  of  resonances  unless  ^ 0-  resonance  occur  for  a value  close  to  (l/c;  in 

such  a case  the  value  of  ap,,  » a Vfl“/c*  - 4^  will  be  very  small  and  the  radiation  insig- 
nificant. unless  the  resonant  frequency  had  been  extremely  large,  (l»nc/a.  The  latter 
is,  however,  impossible  if  the  velocity  of  sound  in  the  shell  and  in  the  medium  are  of  the 
same  order.  The  conclusion  can,  therefore,  be  drawn  that  resonances  can  only  be  sig- 
nificant if  4^—  0. 

Conclusion 

The  essential  meaning  of  Figs.  8a  and  b is  the  fact  that  the  interaction  of  struc- 
tural and  radiation  damping  allows  substantial  radiation  only  in  a relatively  narrow  band 
of  frequenciesi  its  location  and  width  depends  on  the  value  of  n.  The  highest  radiation 
occurs  if  structural  and  radiation  damping  are  about  equal.  Only  if  the  resonant  fre- 
quency of  the  structure  lies  in  this  band  can  noticeable  resonant  effects  ever  occur; 
whether  these  resonant  effects  actually  are  noticeable  depends  in  any  particular  case  on 
the  pressure  level  in  the  other  modes,  and  can  not  be  decided  in  a general  manner.  For 
the  steel  shell  considered  earlier  this  question  is  considered  in  Sec.  VII. 


VI.  Displacements  at  Resonance  and  Effect  of  Rigid  Diaphragms 
a.  Displacements 

Applying  the  reasoning  which  led  to  Eq.  (14)  for  the  pressure  p,  Eq.  (5)  leads 
to  an  expression  for  the  radial  displacement  w due  to  the  loading  Eq.  (16) 

cos  n 9 e***^  f e*^*d^ 

2irmn*  ^ a(ap) 

The  path  of  integration  is  the  same  as  for  the  integral  Eq.  (14),  see  Fig.  2. 

In  general  one  may  expect  that  an  oscillating  load  of  the  type  defined  by  Eq.  (16X 
i.e.  acting  at  z = 0,  will  result  in  displacements  which  are  significant  only  for  a "rea- 
sonable" distance  from  the  force,  comparable  to  the  physical  dimensions  of  the  shell 
(radius  a,  shell  thickness  h).  The  integral  in  Eq.  (32)  was  evaluated  for  resonant  cases 
in  Appendix  6,  in  order  to  show  that  this  is  not  generally  true  for  the  response  of  the 
shell  at  resonance.  The  response  may  then  be  extremely  widespread. 

The  result  of  the  analysis  in  Appendix  B,  Eq.  (B-35),  gives  the  displacement  w 


*From  this  reasoning  it  would  appear  that  the  nearly  complete  suppression  of  radiation 

by  structural  damping  in  similar  circumstances  should  occur  not  only  in  the  case  of 

cylindrical  shells,  but  quite  generally  for  any  type  of  submerged  structure. 


for  resonant  cases  n >■  2 in  the  form  * 


w(z)  = w(o)  e*'**'  (33) 

where  a||  is  in  certain  cases  small  versus  unity.  Consider  a case  * o where 
aflA:  « n;  is  given  by  Eq.  (B-33)  and  due  to  Y„  » 1 one  obtains  (afc,  | ■ n(n  - 1)6 j for 
small  damping  6 and  small  n,  this  gives  a very  small  value  of  . 

The  fact  that  the  amplitudes  of  the  displacement  vary  only  very  slightly  at  dis- 
tances z = a or  even  more,  raises  the  question  of  the  effect  of  a rigid  diaphragm  located 
at  some  distance  z « sa  from  the  load.  In  general,  the  displacement  w.  Fig.  9a,  will  be 
localized,  such  that  a diaphragm  does  not  interfere  with  the  displacements  noticeably; 
no  appreciable  effect  on  the  radiation  can  be  expected.  This  can  no  longer  be  said  if 
the  displacement  is  not  localized,  i.e.  if.  at  resonance,  the  exponent  in  Eq.  (33)  hap- 
pens to  be  small.  One  will  suspect  that  the  diaphragms  will  not  only  prevent  the  dis- 
placement (33),  but  in  doing  so  also  reduce  the  unusually  high  radiation  to  normal  levels. 

The  knowledge  of  the  displacement,  Eq.  (33),  makes  it  possible  to  treat  the  shell 
with  one  or  more  diaphragms  by  applying  the  principle  of  superposition,  and  obtain  the 
effect  of  diaphragms  on  the  radiation  for  the  resonant  cases. 

b.  Effect  of  Diaphragms 

The  effect  of  a " rigid*  diaphragm,  or  stiffening  ring.  Fig.  9b,  at  a distance  sa 
from  the  load  P,  Eq.  (16),  is  essentially  the  application  of  a load  S cos  n 9 e*^at  z = sa 
such  that  the  displacement  at  z = sa  vanishes.**'  The  magnitude,  i.e.  the  absolute  value 
and  the  phase  of  S can  be  determined  from  this  condition.  Observing  that  the  new  force 
is  of  similar  nature  as  P,  Eq.  (33)  gives  also  the  deflections  due  to  S,  if  z is  replaced 
by  z - sa  and  one  obtains 


S = - e‘  (34) 

where 

X = iat,s  (35) 

For  two  symetrically  spaced  diaphragms.  Fig.  9c,  the  forces  at  each  of  them  are 


S 


1 + e 


tx 


(i£i  ) has  a negative  real  part,  which  makes  | s|  < 1 . 


(36) 


*Note  that  w(o)  is  a function  of  t and  0,  it  contains  a factor  cos  n 0 e‘^. 

**This  is  not  true  for  n = 1 where  the  displacement  is  such  that  the  shell  sections  re- 
main circular  and  diaphragms  are  ineffective. 
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Consider  now  the  asymptotic  sound  field  due  to  the  superposition  of  the  original 
force  P and  the  reactions  S due  to  the  diaphragms.  For  points  at  large  distances  r the 
differences  sa  in  location  of  the  forces  P and  S will  be  immaterial,  and  the  entire  effect 
is  that  of  a force  Pj  = P + Z S at  z ■ 0.  For  the  case  of  one  diaphragm. 


and  for  two  diaphragms 


I P,  I - P 1(1  - e*')!-  |k|p 


|P,I  -F 


1 - 


2e^ 

1 + e**^ 


(37) 


where  the  approximations  may  be  used  only  if  |X|  is  sufficiently  small.  If  X is  small  the 
reduction  in  the  force  from  P to  P,  , and  the  corresponding  reduction  of  the  sound  pres- 
sure is  substantial.  The  reduction  for  two  diaphragms,  being  of  the  order  X*  is  much 
larger  than  for  one  only  ; other  arrangements  of  two  or  more  diaphragms.  Fig.  9d,  can 
easily  be  studied;  the  reduction  remains  always  of  the  order  of  X*. 


c.  Reduction  of  the  Resonant  Pressure  by  Diaphragms,  no  Structural  Damping 
For  small  frequencies  n/Z  Eqs.  (26,  30)  give  the  resonant  pressure 


COS  n 0 


irn  V 
~ -r  *n 


(38) 


For  the  case  0 Eq.  (B-33)  gives 


iinjLli 


the  reduction  of  the  force,  Eq.  (37),  reduces  the  pressure  to  | p^  | = ^ jp^^  j,  where  X 
is  given  by  Eq,  (35),  and 


cos  n e 


1=  n(n  - l)s*^ 


(39) 


Yu  being  very  large,  this  is  a considerable  reduction;  comparison  with  Eq.  (27)  shows 
that  the  magnification  n(n  - l)s*  is  only  moderate,  unless  n or  s are  large;  comparison 
with  Eq.  (29)  shows  that  the  reduction  due  to  diai^ragms  for  small  n and  s — 1 is  larger 
than  the  one  due  to  small  structural  damping.  However,  for  larger  n structural  damp- 
ing will  be  more  effective  than  diaphragms. 

For  the  case  » 0,  Eq.  (B-13)  gives  i , being  no  longer  small,  the 

analysis  in  Appendix  B is  invalid.  It  is,  however,  obvious  without  further  analysis  that, 
the  radiation  being  emitted  by  the  mode  of  half  wave  length  L > ir/4,|,  only  diaphragms 
spaced  at  slb-^  L/2  can  possibly  be  effective,  if  a|,l^n  this  requires  s < ir/2n.  Unless 


-16- 


such  close  diaphragms  are  used,  damping  will  in  this  case  be  the  dominant  effect. 

d.  Reduction  of  the  Resonant  Pressure  by  Diaphragms,  including  Structural  Damping 

The  point  of  principal  interest  is  whether  the  peaks  of  the  resonant  pressure 
curves  with  damping,  Fig.  8,  can  be  further  reduced  by  diaphragms  of  reasonable  spac- 
ing, s 1.  Little  further  reduction  below  the  level  with  damping  can  be  expected  if 
» 0,  or  if  n is  large,  even  if  ^ 0,  as  the  preceding  discussion  clearly  shows  that 
damping  in  these  cases  is  relatively  much  more  effective  than  the  diaphragms. 

The  evaluation  shows  that  for  n = 2 the  resonant  pressures  with  damping  are  sub- 
stantially reduced  by  diaphragms:  they  are,  however,  reduced  very  slightly  only,  less 
than  10%,  from  the  pressures  with  diaphragms  but  no  damping.  For  n = 10  the  situation 
is  reversed.  The  pressures  with  damping  and  diaphragms  are  considerably  less  than 
those  with  diaphragms  only,  but  the  reduction  is  only  slight,  about  25%,  if  the  results 
are  compared  with  the  damped  case  without  diaphragms. 

e.  Required  Rigidity  of  Diaphragms 

The  effect  of  " rigid"  diaphragms  or  stiffening  rings,  was  studied  in  this  section 
and  the  question  arises  when  is  a physical  member  stiff  enough  to  be  considered  rigid. 
Common  sense  will  suffice  to  make  the  decision  for  diaphragms,  but  this  is  not  true  for 
stiffening  rings.  The  decision  depends  on  the  relative  values  of  the  frequency  of  the 
ring,  of  the  forcing  frequency  fj,  and  of  the  frequency  of  the  shell  between  the  rings. 

The  situation  is  similar  to  the  one  encountered  in  Sec.  Ill  where  conditions  were  given 
to  ensure  sufficient  stiffness  of  the  shell.  Using  the  symbols  defined  for  Eqs.  (21b,  c), 
stiffening  rings  spaced  at  distances  t = 2sa  may  be  considered  " rigid"  if 

n«  y^l  + M w,(L  = «)  (40a) 

he, 

w,(L  = «)  » 6.2 (40b) 

It  is  interesting  to  note  that  the  frequency  in  vacuo  M|(L  = »)  increases  with  increasing  n, 
and  that  any  stiffening  ring  for  sufficiently  large  n may  be  considered  rigid.  Eqs.  (21b) 
and  (40b)  differ  only  by  the  direction  of  the  Inequality,  and  it  appears  that  for  low  n a 
shell  with  ring  stiffeners  can  be  treated  by  the  ai^roximate  method  of  Sec.  III.  For  large 
n the  same  rings  can  be  considered  "rigid”;  the  structure  can  then  be  studied  first  as  a 
shell  without  stiffeners,  with  a subsequent  correction  for  "rigid"  stiffeners.  While  the 
latter  case  has  not  been  solved  generally,  the  resonant  cases  can  be  treated  as  shown  in 
this  section. 

f.  Conclusion 

The  effect  of  diaphragms  or  rigid  stiffeners  on  the  resonant  pressures  can  only 
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be  large  for  low  values  of  n >■  2,  for  high  values  of  n any  effect  is  insignificant  com- 
pared to  structural  damping.  The  delineation  of  high  and  low  values  of  n depends  on  the 
spacing  of  the  diaphragms.  For  s — 1,  i.e.  a spacing  about  equal  to  the  diameter,  n ■ 10 
is  already  high. 

Even  for  low  values  of  n the  effect  of  diaphragms  is  not  large,  unless  is  small. 
This  is.  however,  of  little  importance  as  in  the  case  of  large  damping  alone  practi- 
cally eliminates  the  effect  of  resonance,  as  was  explained  at  the  end  of  Sec.  V. 

The  effectiveness  of  diaphragms  in  reducing  radiation  at  resonance  will  be  seen 
on  the  numerical  example  in  the  next  section. 


VII.  Pressures  at  Resonance  for  a Typical  Stiffened  Steel  Shell 

The  preceding  two  sections  discussed  the  effect  of  damping  and  diaphragms  in  a 
rather  abstract  manner.  In  this  section  the  actual  resonant  frequehcies  and  the  asso- 
ciated pressure  levels  p will  be  considered  for  a typical  shell  giving  numerical  results 
which  may  be  easier  to  appreciate. 

The  properties  of  the  shell  are;  pa/2m  > 4.  h > aAOO.  stiffener  spacing  i « 0.25a. 
diaphragms  at  distances  of  3a.  making  s > 1,5.  Fig.  10.  structural  damping  6 ■ 0,01. 

The  rigidity  of  the  ring  stiffeners  is  such  that  for  n > 2 the  frequency  in  vacuo  is 
w,(L  * «)  * 0,16  a/c. 

In  the  last  part  of  Sec.  IV  it  was  shown  that  the  analysis  of  the  stiffened  shell  is 
only  applicable  if  aO/c  « 2.7  and  n < 7.  This  range  can  be  studied  with  the  formulae 
developed  so  far.  The  range  can  be  extended  above  these  limits  as  will  be  shown  later 
in  this  section. 

a.  Range  afl/c  « 2.7  and  n < 7 

The  resonant  frequencies  can  be  located  from  Eq.  (22)  which  is  the  frequency 
equation  for  the  submerged  shell  with  A ^ 2/ni  for  (g  ■ 0 this  is  equivalent  to  the  fre- 
quency equation  for  a ring  carrying  in  addition  to  the  mass  m of  the  shell  also  the  vir- 
tual mass  of  fluid  (see  Ref.  1.  p.  12) 

my-^«,A — (41) 
n*+  1 

where  A — 2/n  and  a|—  n*/n*  1.  The  frequencies  of  a ring  in  vacuo  vary  as 

, 

and  as  the  inverse  of  the  square  root  of  the  mass.  Tying  into  the  value  0.16  a/c  in  vacuo 
for  n • 2 gives  the  resonant  frequencies 
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an  , 0.060  n(n*  - l)a  ^^2) 

c|/n«+l+^ 

V in 

Numerical  values  of  n are  listed  in  Table  I.  It  is  noted  that  Eq.  (22)  in  this  case  has  no 
root  for  n > 1 ; there  are  therefore  no  resonant  frequencies  for  n = 0 or  1. 

Using  Eqs.  (26,  28.  35,  37,  and  B-33)  the  resonant  pressures  with  and  without 
damping  and  diaphragms  were  computed  (with  slide  rule  accuracy)  for  the  case  > 0, 
i.e.  at  z > 0,  and  are  shown  in  Table  I.  All  frequencies  are  sufficiently  low  so  that  the 
approximation  Eq.  (B-30)  could  be  substituted  into  Eqs.  (26,  28). 


Table  I.  Pressure) p|  at  Resonance  at  z ■ 0 


«r  IpI  /cos  n 0 * 

n 

Resonant 

Frequency 

an/c 

-Y« 

no  damping 
no 

diaphragms 

with  damping 
no 

diaphragms 

no  damping 
with 

diaphragms 

with  damping 
and 

diaphragms 

2 

0.08 

200 

630 

0.50 

0.022 

0.022 

3 

0.25 

330 

780 

0.30 

0.041 

0.041 

4 

0.52 

426 

2680 

0.24 

0.063 

0.063 

5 

0.9 

0.25 

0.103 

0.103 

6 

1.4 

0.28 

0.19 

0.19 

The  diaphragms  are  clearly  dominant  in  determining  the  level  of  the  pressure  in  this 
range  as  the  figures  with  diaphragms  and  damping  are  indistinguishable  from  those  with 
diaphragms  but  no  damping. 

To  judge  the  importance  of  the  pressure  p reduced  by  the  diaphragms  it  should 
be  compared  with  the  pressures  p for  the  other  values  of  n for  the  same  frequency. 

Fig.  6a.  It  must  be  kept  in  mind  that  only  the  total  pressure  comprising  the  sum  of  the 
components  for  all  values  of  n can  be  observed,  and  that  a peak,  say  in  the  n > 6 com- 
ponent, will  be  lost  in  the  total  pressure,  if  this  peak  was  small  in  comparison  to  the 
pressure  in  the  modes  n > 0, 1 and  2.  This  is  the  position  in  this  range  for  all  n in  case 
of  the  shell  with  diaphragms;  it  can  be  predicted  that  regardless  of  the  exact  nature  of 
the  applied  force  the  total  pressure  will  not  show  appreciable  signs  of  resonant  frequen- 
cies. (This  statement  excludes  the  unlikely  case  of  a force  which  contains  only  the  n^ 
component.) 

If  the  shell  has  no  diaphragms  the  position  is  different;  the  case  n = 2 would  pro- 
duce a qtiite  pronounced  resonance,  as  the  level  of  the  pressure  is  about  10  times  larger 
than  the  largest  of  the  other  components.  The  result  becomes  then  very  sensitive  to 
changes  in  the  value  of  structural  damping. 

Table  I also  shows  that  with  increasing  n the  resonant  pressure  levels  increase 
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rapidly.  This  is  in  part  due  to  the  fact  that  the  diaphragms  become  less  effective  for 
larger  n.  The  damping  effect  is  also  peculiar,  as  the  pressure  level  decreases  from 
n > 2 to  5,  but  then  increases  again;  this  is  due  to  the  shape  of  the  curves  in  Figs.  8. 

The  peak  of  the  damped  pressure  curve  moves  gradually  towards  increasing  values  of 
n,  while  the  resonant  frequency  increases  much  faster,  nearly  like  n*.  The  resonant 
frequency  therefore  gradually  moves  into  the  peak  of  the  response  curve,  which  accounts 
for  the  increase  of  | p | for  n > 5.  Important  resonances  might  occur  for  the  shell  under 
consideration  above  n ■ 6 requiring  further  consideration. 

b.  Refined  Theory  for  Higher  Values  of  fl  and  n 

The  three  limitations  Eq.  (21a  - c)  originate  from  the  basic  assumption  that  the 
deflection  of  the  shell  between  stiffening  rings  must  be  negligible  compared  to  the  dis  - 
placements  of  the  rings.  It  is  possible  to  approach  the  problem  without  this  assumption 
in  the  following  manner. 

Let  the  deflection  of  the  shell  between  stiffening  rings  be  proportional  to 
1 - cos  2vz// ; the  entire  deflection  is  then 

w ■ cos  n 8 e***^sin  |z  + B^l  - cos 

where  B is  an  arbitrary  value.  Provided  the  distance  between  nodes  in  the  longitudinal 
direction,  t/Z,  is  smaller  than  the  distance  between  circumferential  nodes,  ira/x  the 
shell  in  vacuo  will  vibrate  between  stiffening  rings  with  a frequency  roughly  that  of  a 
clamped  plate  of  span  t 

he, 

„,~6.2_  (44) 


The  above  condition  on  the  nodes  requires 


n « 


2wa 

/ 


(45) 


For  the  above  example  this  gives  n « 8ir,  a considerable  extension  of  the  range. 

Use  is  now  made  of  the  fact  that  each  resonant  frequency  n is  essentially  the  real 
part  of  a complex  frequency  of  free  vibration  of  the  submerged  shell;  the  resonance  be- 
ing pronounced  only  if  the  imaginary  part  of  the  frequency  is  small.  The  real  parts  of 
slightly  complex  frequencies  of  free  vibrations  can  be  found  as  frequencies  of  the  shell 
in  vacuo  if  the  actual  mass  is  increased  by  an  appropriate  virtual  mass  of  the  medium 
(Ref.  1,  p.  10-12).  The  problem  of  finding  the  resonant  frequencies  is  therefore  reduced 
to  the  one  of  finding  frequencies  of  a stiffened  shell  in  vacuo  having  a slightly  different 
mass  (This  approach  was  already  used  in  deriving  Eq.  (42) ). 

Having  made  the  assumption  (43)  concerning  the  deflections  of  the  shell  between 


-20- 


the  stiffeners,  the  combined  system  of  shell  and  rings  can  be  treated  by  the  method  of 
Ref.  5.  The  resonant  frequencies  n,  including  the  effect  of  the  shell  deflection,  are  the 
roots  of  the  equation 


(46) 


where  n is  the  resonant  frequency  if  the  shell  is  rigid,  Eq.  42 ; (Ig  is  the  frequency  of 
the  shell  vibrating  between  stiffeners.  Eq.  44,  but  with  a correction  for  the  effect  of  the 
virtual  mass 


Os 


(47) 


m„  is  the  mass  of  the  stiffeners  per  unit  length;  m^is  given  approximately  by  Eq.  (41). 

To  defend  the  procedure  it  is  necessary  to  consider  the  interaction  problem  be- 
tween shell  and  medium,  treated  in  Ref.  1 for  the  case  without  local  shell  deformation. 
The  entire  displacement  Eq.  (43)  may  be  written  in  the  form 

w ■ cos  n 9 e*®*  sin  + B,  sin^  + "»•  B,  sin  (48) 

where  W,  B|  and  B|  are  new  constants.  This  equation  being  a sum  of  sine  terms  in  i the 
fundamental  procedure  of  Ref.  1 can  be  applied,  but  there  will  be  coupling  between  the 
modes  having  the  parameters  ( and  ( ± 2v//.  While  this  complicates  matters  it  can  be 
expected  not  to  change  the  results  qualitatively. 

As  far  as  resonant  radiation  is  concerned  it  should  be  noted  that  if  ( ■ fin  ~ 0 and 
ji  is  small  the  values  i ± iTt/i  will  be  so  large  that  these  modes  do  not  radiate  at  all,  and 
the  intensity  of  radiation  will  be  given  by  the  various  formulae  derived  earlier.  This 
can  be  expected  if  the  condition 


an  . aw 

c ^ / 


(49) 


is  satisfied  in  which  case  the  suggested  procedure  will  give  particularly  good  results. 

It  is  therefore  believed  that  the  results  concerning  resonant  pressures  obtained 
in  Secs.  V and  VI  remain  qxialltatively  applicable  if  the  resonant  frequencies  are  ob- 
tained from  Eq.  (46)  provided  Eqs.  (4S)  and  preferably  also  (49)  are  satisfied. 

The  two  resonant  frequencies  n,  ^of  the  shell  considered  were  obtained  from 
Eq.  (46)  and  are  shown  in  Table  II;  for  comparison  the  frequency  n neglecting  the  effect 
of  shell  deflection  between  rings  is  also  shown. 
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Table  II.  Resonant  Frequencie»  of  Stiffened  Steel  Shell 


n 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

an/c 

0.08 

0.25 

0.52 

0.90 

1.40 

1.95 

2.65 

3.5 

4.4 

5.4 

6.6 

7.9 

aii,  /c 

0.08 

0.25 

0.51 

0.86 

1.23 

1.61 

1.86 

2.10 

2.28 

2.38 

2.50 

2.65 

an,/c 

2.7 

3.0 

3.3 

3.7 

4.1 

4.6 

5.3 

6.5 

7.1 

9.3 

11.1 

13.0 

For  low  n the  revised  frequency  equals  the  original  one,  Hi  the  second  frequency 
which  did  not  appear  in  Table  I is  so  high  that  it  lies  outside  the  range  of  pronounced 
resonance  a(l/c4  n/Z.  For  n > 6,  the  lower  values  H,  lie  in  the  range  where  damping 
practically  prevents  radiation,  as  can  be  seen  from  Fig.  Bb  for  n ■ 10.  The  higher  fre> 
quencies  n,.  on  the  other  hand,  are  important  because  they  lie  for  n 8*10  in  the  peak 
zone  of  the  pressure  curve  with  damping.  Fig.  8b.  For  larger  values  of  n the  value 
aili/c  is  larger  than  n and  no  pronounced  resonance  is  possible. 

In  the  range  n ~ 10  diaphragms  reduce  the  resonant  pressure  only  sli^tly  and 
relatively  high  resonant  pressures  in  the  components  n ■ 8 to  10  must  be  expected  at 
their  respective  resonant  frequencies.  For  6 • 0.01,  which  is  quite  a low  damping  value, 
the  level  would  be  about  5 times  that  of  the  other  components  for  smaller  n.  However, 
if  one  takes  into  account  that  there  are  many  of  the  other  components  which  will  add  up^ 
and  that  the  expansion  coefficient  of  the  force  for  high  n is  usually  lower  than  for  low 
n,  one  can  not  expect  pronounced  signs  of  these  resonances  in  the  total  pressure.  There 
might  be  a doubling  of  the  pressure,  but  no  hi|^  peak.  If  the  damping  is  higher  than  the 
assumed  6 > 0.01,  the  effect  would  even  be  less. 

The  refined  analysis  just  discussed  is  subject  to  the  limitations  Eqs.  (4S)  i.e.  in 
the  presently  discussed  case  n < 8ir;  for  the  range  of  significant  results  n ^ 8-10,  limi- 
tation (49)  giving  afhb  < 4v  is  also  satisfied.  This  indicates  that  resonances  for  n < 8ir 
have  essentially  been  covered. 

p.  The  Range  n » 8w 

For  high  values  of  n the  approach  mentioned  in  Sec.  Vie  can  be  used;  the  stif- 
fening rings  are  for  high  values  n so  rigid  that  they  may  be  treated  as  diaphragms.  The 
resonant  frequencies  of  the  actual  stiffened  structure  are  then  those  of  the  shell  without 
any  stiffeners,  only  the  magnitude  of  the  resonant  pressures  being  affected  by  the  exist- 
ence of  the  'rigid*  stiffeners. 

Considering  the  case  ~ 0,  the  resonant  frequencies  of  the  shell  are  those  of  a 
ring  of  thickness  h;  for  large  n the  virtual  mass  of  the  medium  becomes  negligible  and 
the  frequency  will  be  approximately  equal  to  that  of  a flat  plate  having  a span  J equal  to 
the  distance  between  the  nodal  generators.  T • aVn. 

he. 


The  resonant  pressure  | | allowing  for  damping  for  large  values  of  n is  a func  - 

tion  of  the  frequency  similar  to  the  one  shown  in  Fig.  8b  for  n » 10,  except  that  the  peak 
will  be  closer  to  n and  narrower.  This  means  only  those  resonances  will  not  be  elimi- 
nated by  damping  for  which  the  value  an/c  given  by  Eq.  (50)  is  slightly  less  than  n.  The 
interesting  values  of  n are  therefore 


n-3.5^  (51) 

For  the  case  under  discussion  a/h  » lOO.  c,  /c  ^ 3.5,  giving  n » 100  and  a(l/c  ~ 100. 
Due  to  the  high  value  of  n.  the  distance  between  nodal  lines,  avAOO,  is  so  much  smaller 
than  the  distance  between  the  ring  stiffeners  A > a/4,  that  the  rings  can  have  no  appre- 
ciable effect.  The  expected  resonant  pressure  can.  therefore,  be  estimated  by  Eq.  (30), 

The  pressures  p of  the  other  components  have  a level  given  by  Eq.  (31),  for  the  selected 
shell  p ~ 1.  As  there  are  nearly  a hundred  of  these,  the  total  pressure  will  again  not 
show  a pronounced  sign  of  the  resonance. 

d.  Conclusion 

Except  for  transition  points  the  entire  range  has  been  covered,  but  no  indication 
of  pronounced  resonance  in  the  pressures  at  large  distances  were  found  for  the  steel 
shell  in  water  selected  as  example.  Mild  resonant  effects  increasing  the  total  pressure 
(to  maybe  twice  its  regular  value)  can  be  expected  in  the  range  afl/c  ^ 6 and  again 
an/c  ~ 100. 

The  analysis  did  not  indicate  that  the  results  are  sensitive  to  changes  in  the  di- 
mensions of  the  shell  i no  pronounced  resonant  effects  should  therefore  be  expected  for 
any  shell  of  this  type. 


Vin.  Resonances  for  n » 1 

Eq.  (22),  which  gives  the  location  of  resonant  frequencies,  has  in  this  case  no 
solution  at  all  for  the  steel  shell  in  water  considered  as  example.  One  can  find  such 
solutions  for  steel  shells  with  thicker  walls,  but  only  for  much  heavier  shells,  if 
pa/Sm4  )/2.  However,  even  these  solutions  pertain  to  values  which  are  nearly  equal 
to  n/c.  i.e.  to  cases  which  are  for  practical  purposes  eliminated  by  damping,  as  ex- 
plained at  the  end  of  Sec.  V. 

In  the  case  n • 1 no  solution  of  Eq.  (22)  exists  ever  for  steel  shells  at  or  near 
» 0.  while  for  n ^2  solutions  could  be  found  near  (r  * Oj  these  latter  solutions 
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caused  the  widespread  and  relatively  pronounced  resonances  for  the  shells  without  dia- 
phragms. The  different  behavior  for  n ■ 1 is  caused  by  the  fact  that  for  n ^2  the  sub- 
merged shell  has  a slightly  complex,  low  frequency  of  free  vibration  associated  with  a 
mode  in  which  the  generators  of  the  cylindrical  shell  remain  straight i the  correspond- 
ing frequency  for  n > 1 degenerates  into  (1  = 0,  a rigid  body  motion  of  the  shell,  and  the 
resonance  disappears. 

If  one  introduces  in  some  way  a restoring  force  which  would  resist  rigid  body 
motions  of  the  shell,  the  trivial  frequency  (1  ■ 0 would  disappear  and  resonances  in  the 
mode  n > 1 may  occur.  There  is  an  interesting  possibility  of  introducing  such  a restor- 
ing force.  Consider  a floating  rigid  shell  with  a mass  center  below  the  center  line  of 
the  cylinder  in  a gravitational  field ; the  shell  can  execute  oscillations  (rolling)  during 
which  the  circular  cross-section.  Fig.  11,  will  be  displaced  laterally,  coupling  the  body 
motion  with  the  n « 1 oscillations  of  the  fluid.  Such  a system  will  have  resonances  of 
the  type  found  previously  for  n ^ 2.  This  is  presumably  also  true  if,  instead  of  the  rigid 
cylinder,  an  elastic  shell  with  eccentric  mass  is  considered.  It  may  well  be  that  the 
resonant  radiation  in  this  case  would  be  quite  high,  as  diaphragms  will  have  no  effect;* 
furthermore,  the  rigid  body  part  of  the  motion  of  the  shell  will  have  no  structural  damp- 
ing. This  is  replaced  by  the  viscous  damping  due  to  the  relative  motion  of  shell  and 
medium,  but  may  be  very  low. 


IX.  Qualitative  Validity  of  Results  for  Lx>ng  Shells  of  Finite  Length 

If  the  displacements  of  an  infinitely  long  shell  due  to  an  exciting  force  are  local- 
ized as  shown  in  Fig.  9a,  one  can  reason  that  the  sound  field  can  not  be  affected  mate- 
rially if  parts  of  the  shell  at  large  distances  from  the  force  are  replaced  by  the  acoustic 
medium.  The  results  found  for  the  infinite  shell  can,  therefore,  be  applied  to  long  but 
finite  shells,  if  the  load  does  not  act  close  to  the  end,  except  in  those  cases  where  the 
displacements  are  not  localized.  This  excludes  certain  resonant  cases,  and  also  the  re- 
sponse in  the  n > 1 mode  for  low  frequencies.** 

If  the  shell  has  a number  of  rigid  dia{4tragms  one  can  reason  further  that  their 


*Sm  footnote  on  page  15. 

**  The  displacement  of  an  infinitely  long  bar  in  vacuo  due  to  a transverse  oscillating 
force  is  w = e'^*^  where 


K - 


v: 


El 


a*mn* 


For  sufficiently  small  frequency  (1  the  decay  constant  X will  become  small,  and  the 
displacement  can  no  longer  be  considered  localised.  (E  is  Young's  modulus^  I the 
moment  of  inertia  of  the  bar.) 
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existence  enforces  localized  displacements,  excepting  the  n ■ 1 node;  in  particular  the 
conclusion  can  be  drawn  that  such  shells  do  not  have  pronounced  resonances  for  n ^ 1. 

No  such  conclusion  can  be  drawn  for  n ■ 1,  as  the  infinitely  long  shell  for  low 
frequencies  can  not  be  used  as  approximation  for  a shell  of  finite  length. 
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APPENDIX  A 


Evaluation  of  the  Integral  In  Eg.  (14)  for  Large  Values  of  r 

The  integral  in  Eq.  (14)  may  be  conaidered  in  the  complex  plane  and  can  be 
evaluated  along  a suitable  selected  alternative  path.  When  distorting  the  path  one  must 
consider  that  the  restrictions  in  Eq.  (9)  on  p impose  a restriction  on(,  and  it  is  easily 
seen  that  ( = | 4 | is  only  permitted  to  assume  values  such  that  0 < 0 < ir/2  or 
ir  ^ 9 4 3w/2  . 

The  original  path  of  integration  as  shown  in  Fig.  2 is  interrupted  at  each  pole  . 
This  path  is  altered  in  two  ways . Fig.  A -1 ; 

1 . It  is  made  continuous  by  adding  semi  -circular  paths  of  radius  c above  the 
real  axis  for  positive  poles  and  below  for  negative  ones.  This  addition  in- 
creases the  value  of  the  integral  by 


where  Rj  is  the  residue  of  the  integral  at  , 

2.  The  path  near  the  branch  points  is  distorted  to  a semi -circle  of  radius  7. 
As  the  path  remains  on  the  proper  side  of  the  branch  cut  the  value  of  the 
integral  does  not  change;  if  a pole  occurs  within  this  circle  its  residue  is 
already  accounted  under  1. 

The  value  J of  the  integral  in  Eq.  (14)  is 


J 


f H‘*’(pr)G(t)e‘«*  di  „ 

J „<*»  Zm 

-iNAla^J  j 


(A-1) 


where  the  sign  under  the  sum  is  the  same  as  that  of  the  respective  value  • 

The  use  of  a complex  path  introduces  a difficulty  concerning  the  function  N de- 
fined by  Eq.  (6).  This  equation,  taken  from  Ref.  1,  applies  only  for  real  values  of(. 
because  it  refers  to  the  frequencies  W|,  of  modes  of  the  cylindrical  shell  in  vacuo  having 
a half  wave  length  L » I v/t  I . If  i becomes  complex  this  definition  loses  its  meaning. 
However,  one  can  visualise  a solution  of  the  problem  under  consideration  by  integral 
transform,  in  which  case  the  solution  would  contain  an  integral  similar  to  Eq.  (A-1). 
Such  a solution  would  give  a more  general  definition  of  the  function  N valid  also  for 
complex  values  of  It  is  therefore  assumed  that  N may  be  considered  as  a function  of 
the  complex  variable  (,  althou^  only  the  values  for  real  ( can  be  computed  from  Eq. 
(16).  It  will  be  seen  later  that  values  of  N for  complex  argument  will  not  be  required  for 
the  present  purposes. 

Consider  the  values  of 
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along  the  new  path.  The  smallest  values  of  |P|  occur  when  4 goes  through  the  semi- 
circles in  the  vicinity  of  4 * ± Ok.  For  small  values  of  7 the  absolute  value  | P|  has  a 
lower  bound 


min  I p I > 


(A -2) 


This  being  the  case,  one  can  select  a radius  R large  enough  so  that  for  all  values  of 
r > R the  asymptotic  expression 


Urn  hI,**  (Pr)  « e * ® 

r -•« 


Vt^ 


(A-3) 


may  be  used.  Substituting  into  Eq.  (A-1) 


l‘mJ  = ,yiJ*Mr  (j(U-fr) G(4)d4 

*'^*  ''*  J . r. 


(pa)  [l  - 


Z*-  " r 


(A -4) 


This  integral  is  of  the  form 


I . f e^^‘*  G(4)d4 

•^A 


where  4 = x -t-  iy  is  a complex  variable,  and 


G(4) 


Giil 


^Hi‘*(pa)  [l  - jsifsjipj] 

f(4)  = iy4  - ip  = i 


(A -5) 


(A -6) 


(A-7) 


(A -8) 


Integrals  of  the  type  (A -5)  can  be  evaluated  approximately  for  large  values  of  r 
by  utilizing  the  fact  that  the  significant  contribution  to  the  value  of  the  integral  is  due  to 
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the  portion  of  the  path  in  the  vicinity  of  "saddle  points"  40(that  is  points  where  f'(4«)  " 0), 
provided  G(()  has  no  singularity  on  the  path. 

The  condition 


f'(6o)  “ i 


(A -9) 


furnishes  two  roots 


^0 


/l  + y* 


of  which  only  one  pertains  to  the  branch  of  the  function  t{i)  applicable.  Eq.  (A -9)  may 
be  written 


Y + ^=0  (A-9a) 

and  because  of  the  restrictions  in  Eq.  (9),  p must  be  real  and  not  negative  if  | is  real. 
Assuming  z > 0.  y > 0 Eq.  (A -9a)  excludes  positive  values  of  and  the  only  saddle 
point  is  at 


= *0+ iyo-  *0*  - ?-F=^^==-*  yo  ' 0 (A -10) 

The  selected  path  goes  through  this  saddle  point  and  follows  in  its  vicinity  the 
real  axis.  The  function  f(4)  > 4 i^  is  purely  imaginary  for  real  values  ( ^ x of  the 
argument.  The  path  being  such  that  the  real  part  4 " 9,  the  situation  is  suitable  for  the 
application  of  the  method  of  stationary  phase,  (Ref.  3,  p.  505,  506).  For  large  values  of 
r the  integral  I can  be  approximated  (Ref.  3,  Eq.  (9)) 

m/  r -|  G(i.)exp[inKe,)+^l  +oQ  (A-11) 

where  G and  4e  defined  by  Eqs.  (A -6)  and  (A-10).  while 

*(«,)»  Vi0-Ps  (A -12) 

if'iU  « (A-13) 


Po* 


(A -14) 
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The  square  root  in  Eq.  (A -14)  must  be  taken  as  positive. 
Eq.  (A -4)  becomes 


lim  J _ 
r -•  CO 


2G{e,)  exp  i(lli!^  - «*  )] 

Vr*  + z*  Hi*'  (aPo)  [l  - 


lim  Rj 

r -CO 


(A -15) 


where  the  term  containing  0(l/r)  is  omitted.  The  residual  terms  are  still  to  be  deter- 
mined. 

The  term  H*||*'(pr)  in  Eq.  (14)  governs  the  behavior  of  the  residues  Rj . As  was 
stated  previously  poles  occur  only  for  1 4j  | > fl/c . We  consider  first  the  case 

(A -16) 

as  the  case|4j|  = fl/c  requires  special  treatment.  Eq.  (0)  indicates  that  P = - (iy)  is  a 
negative  imaginary  number;  the  function  h|,*'  (-iyr)  decays  as  /Yt , which  is  faster 
than  the  first  term  of  Eq.  (A -15),  the  asymptotic  behavior  of  which  depends  on  the  term 

/r^+z* 

in  the  denominator.  Residues  for  poles  |ij  | > (2/c  give  therefore  no  contribution. 

If  lij  I n/c  the  coefficient  P - 0,  and  the  behavior  of  the  residues  is  governed  by 
the  value  of 


11„  h'"  (w 

? - 0 h|,"  (pa) 


(A -17) 


If  n ^ 2 the  residues  decay  faster  than  the  first  term  of  Eq.  (A -15),  but  the  cases  n » 0 
and  n = 1 still  remain. 

Even  the  cases  n = 0 and  1 do  not  give  contributions  because  the  values  Rj  of  the 
residues  vanish  if  the  poles  are  at  ± nk.  This  fact  can  be  shown  by  a routine  computa- 
tion not  presented  here,  and  is  due  to  the  behavior  of  Ho"  and  h{*'  and  their  derivatives 
for  vanishing  argument.  The  residues  never  give  a significant  contribution  and  the  first 
term  of  Eq.  (A -15)  represents  the  entire  result. 

To  obtain  an  estimate  concerning  the  values  of  r for  which  the  result  is  appli- 
cable one  must  consider  the  two  approximations  made.  The  principal  approximation 
lies  in  Eq.  (A  -11)  which  indicates  the  second  neglected  term  is  of  the  order  |/a7r  smaller 
than  the  first.  This  means  r/a  > 50  will  give  some  measure  of  approximation.  The  sec- 
ond a{^roximation  is  that  Eq.  (A -3)  is  usable*  only  if  pr  > 2n  1.  Enlarging  the  radius 


*Eq.  (A-3)  gives  for  n<  6 less  than  5%  error  in  the  absolute  values  of  H(z)  if 


c of  the  semi -circle  around  the  branch  point  in  Fig.  A-1.  one  can  go  only  as  far  as  the 
saddle  point,  and  min  p according  to  Eq.  (A -2)  must  be  less  than 


This  gives  the  condition 


rPo  = — - — > 2n  + 1 (n  < 6)  (A-18) 

Comparing  the  two  conditions  for  z = 0 one  finds  Eq.  (A-18)  can  usually  be  disregarded, 
the  exception  being  quite  low  frequencies  fi.  Eq.  (A -18)  controls  only  if  an/c  0.02  (2n  + l). 


(footnote  continued) 

I z I > 2n  -«■  1;  the  error  in  the  phase  may  be  larger  but  is  of  no  consequence  as  it  doe^ 
not  affect  the  stationary  phase  method  as  such,  and  we  are  ultimately  interested  in  | p|. 


APPENDIX  B 


Estimate  of  the  Radial  Diaplacement  w at  Reaonance,  n » 2 

The  purpose  of  this  appendix  is  to  demonstrate  that  at  resonance  the  displace- 
ments w may  be  substantial  at  much  larger  distances  from  the  applied  load  than  one 
would  normally  expect.  It  is  intended  to  apply  the  results  to  stiffened  shells;  for  such 
shells  the  function  N is  nearly  constant  * for  | ( | < n/a.  The  following  computations  are 
restricted  to  the  case  N > constant  for  1 4 | < n/a. 

The  information  desired  is  of  qualitative  nature  and  rather  crude  approximations 
were  considered  permissible. 

The  integral  in  Eq.  (32)  to  be  determined  is 


There  are  no  poles  on  the  path  of  integration;  the  integrand  has  branch  points  at  + fl/c 
and  the  path  passes  below  the  negative  but  above  the  positive  one. 

The  plan  of  action  is  as  follows:  in  cases  of  resonance  the  integrand  1 becomes 
very  large  in  the  vicinity  of  4 = ± 4*  where  4n  defines  the  directions  of  the  largest  radi- 
ation, see  Secs.  V and  VI.  These  large  values  are  caused  by  poles  in  the  vicinity  of  + 4r 
Deducting  a function 


f(e) 


e-  e. 


from  the  integrand,  Eq.  (B-1)  becomes 


(B-2) 


J » 


f)d4  + 


(B-3) 


where  Rj  are  the  values  of  the  residues.  The  integrand  of  the  first  integral  will  nowhere 
be  large,  and  the  integral  will  be  small  compared  to  the  second  one.  The  integrand  of 
Eq.  (B-1)  being  a symmetric  function  of  x.one  can  restrict  the  computation  to  s > 0,  and 
the  second  integral  can  be  evaluated  in  an  obvious  manner  by  the  residue  method. 


J(|z|)  = 2 


[I  - f]d4 


(B-4) 


*lf  the  shell  is  not  stiffened  the  approximation  N > constant  can  also  be  made,  but  is 
valid  only  for  much  lower  values  of  1 4 I . 
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where  only  reeiduei  for  Imtj  > 0 ere  to  be  included. 

Determination  of  the  Poles  and  Reelduee 

The  poles  are  in  the  vicinity  of  |4n|  < n/a  and  it  will  be  assumed  that  N in  this 
range  can  be  treated  as  a constant.  Expressing  the  function  A in  terms  of  J and  Y func 
tions  (as  in  Sec.  V)  the  integrand  becomes 


where 


(B-5) 


The  Bessel  functions  J and  N have  the  argument  aPn  and  the  primes  indicate  deriva- 
tives with  respect  to  the  argument. 


The  subscripts  n will  be  omitted  in  the  following  for  simplicity. 

At  resonance  the  otherwise  dominant  imaginary  part  of  the  denominator  of  Eq. 
(B-5)  vanishes  for 


i 


±6n.P  - Pa 


this  gives  the  relation 


- m ' 0 (B-6) 

where  the  notation  Y(R)  sUnds  for  Y(aPn). 

It  is  inherent  that  in  the  resonant  range  Y and  Y*  » J and  J'  j J and  J'  will,  there 
fore,  be  neglected  versus  Y and  Y',  but  not  versus  products  AY  and  AY'.  Drooling  the 
last  term  in  Eq.  (B-6),  this  equation  is  used  to  eliminate  N in  Eq.  (B-5), 

^ . m UPY'Y'(R|^- 

where 


« J Y'(R)  - Y(R)  J'  ^ - AY(R)  Y'  + i J-  YY*(R)  + Y(R)  Y'  ^ (B-8) 
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Negligible  terms  which  will  not  be  carried  further  have  been  crossed  out. 

To  find  the  poles,  let  D ■ D(a^)  « 0.  Knowing  that  the  roots  occur  near  aPn  the 
function  D is  expanded  in  a power  aeries  in  terms  of 

« « ap  - ap,  (B-8a) 

giving 

D - D(R) +«D'(R)  + . , . » 0 (B-9) 

where  D(R)  > D(aPn).  Using  the  first  two  terms  of  the  series  only,  and  noting  the  rela- 
tion for  the  Wronskian.  (see  Eq.  (25)) 

D » JY'  - YJ'  - 4YY*  - - 6YY'  (B-10) 

*»Ps 


where  all  functions  are  of  argument  apn  (the  notation  J(R)  is  no  longer  necessary).  In 
determining  D'  it  is  found  that  the  entire  contribution  from  the  real  terms  is  small 
versus  the  one  from  the  imaginary  ones,  and  only 

D'  = i [-  Y'*  YY<  + YY"j  * i [y‘* 

remains.  The  root  « of  Eq.  (B-9)  is 


€ 


(B-12) 


It  is  necessary  to  point  out  that  the  denominator  of  Eq.  (B-12)  changes  sign  as  the  argu- 
ment aPn  goes  from  zero  to  n,  and  will  be  0 somewhere;  in  the  vicinity  of  this  point  the 
neglected  terms  will  dominate.  It  can  be  shown  that  this  occurs  at  a point  where  J is 
no  longer  very  much  smaller  than  Y,  that  is,  it  occurs  outside  the  resonant  range. 

Solving  Eq.(B-8)  for  the  values  of  i corresponding  to  the  root  < gives  the  location 
of  the  poles 

D-Ui  (B-U) 


Only  the  pole  i , which  lies  near  has  a positive  imaginary  part  and  is  to  be  included 
in  Eq.  (B-4).  The  required  residue  R|  of  the  integrand  I,  Eq.  (B-7).  at  ( " is 
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Noting  that  « is  small  the  argument  ap  of  tiie  various  functions  can  approximately  be 
taken  as  ap„;  D'  is  then  given  by  Eq.  (B-U)t  using 


one  obtains 

R,  ■ Ae‘^** 


where 
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Estimate  of  the  Remainder  Integral  in  Eq.  (B-4) 

It  is  assumed  that  the  damping  term  in  Eq.  (B-1)  can  be  neglected.  Asymptotic 
values  for  large  £ will  be  required  for  f(4)  and  N.  Noting  and  R^  = -Ae*'^<‘  one 

obtains  from  Eqs.  (B-2,  IS.  16) 


2fi,cos(,z  2i  sin  z 


(B-17) 


The  asymptotic  value  of  N is  obtained  from  the  following  consideration;  for 
4 the  radial,  tangential  and  longitudinal  vibrations  (or  waves)  in  a thin  shell  will  be 
essentially  like  plane  waves  in  a flat  plate,  and  will  not  be  coupled.  The  values  tf|^, 

Eq.  (7){  depend  on  the  radial  component  of  the  vibration^  due  to  the  uncoupling  <T|  = 1 for 
the  radial  (i.e.  bending)  vibrations,  and  ■ 0 for  the  others.  The  bending  frequency 

can  be  written  w,  « c,4  where  c,  is  the  velocity  of  propogation  of  these  waves,  i.e.  the 
velocity  of  shear  waves  in  the  shell  material.  Eq.(B-6)  gives 
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where  fl*  was  neglected  versus  w*  . 

In  the  beginning  of  the  appendix  it  was  assumed  that  N is  constant  for  small  4; 
the  value  of  N is  implied  by  Eq.  (B-6) 
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the  approximation  being  permissible  as  the  argument  of  the  Bessel  function  is  less  than 
n,  (see  Eq.  (B-30) ).  There  is  necessarily  a transition  between  the  extreme  values  given 
by  Eq.  (B-18,  19);  however,  as  an  approximation  it  will  be  assumed  here  that  these  rela- 
tions are  valid  up  to  the  intersection  of  the  two  curves  at  a value  4 * ± 4|, . 

6t  ■ ^ (B-20) 


The  remainder  integral  may  be  written,  noting  1(4)  ■ I(-4) 


(B-21) 


It  is  now  reasoned  that  near  the  origin,  i.e.|  4|  < the  integrand  (I  - f)  must  be  of  the 
order  of  magnitude  of  N,  as  the  extremely  high  values  which  occur ed  due  to  the  cancel- 
lation of  the  a and  N terms,  have  been  removed.  The  third  integral  is  therefore,  of  the 
order 


The  first  and  second  integral  can  be  found  by  substituting  Eq.  (B-17) 


(B-22) 


St. 


The  fourth  integral  will  have  its  largest  value  for  z > 0,  giving  the  limit 


(B-23) 


4 £ 

Collecting  the  variotis  terms  the  value  of  J becomes 


(B-24) 


J(z)  ~ 2wiA  e'*'* 


(B-25) 


valid  if  the  first  term  is  dominant,  and  provided  that  the  terms  originating  from  Eqs. 
(B-22,  24)  are  small 
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It  is  Intended  to  apply  this  result  only  in  cases  where  is  very  smalli  both  terms  of 
Eq.  (B-25)  represent  in  such  cases  functions  which  are  nearly  constant  for  | z|^  a,  and 
it  is  this  fact  which  is  the  point  at  issue.  For  the  present  purpose  the  second  term  of 
Eq.  (B-25)  which  does  not  change  the  character  of  the  solution  can  therefore  be  neglected, 
even  if  it  is  much  larger  than  the  terms  excluded  by  Eq.  (B-26).  The  value  of  J becomes 
finally 

J(z)  --  2wiA  e ' (B-27) 


where  the  condition  (B-26)  is  vital,  while  it  is  sufficient  to  require  less  stringently 
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Simplifications  for  Low  Frequencies 

Useful  simplifications  are  possible  if  the  frequency  is  low  aO/c n/2 ; these 
approximations  cover  the  most  interesting  range  of  high  radiation,  but  are  not  suitable 
for  the  transition  range  of  gradually  diminishing  radiation  as  aO/c  increases  toward  n. 
The  dominant  term  in  the  expansions  of  Yn  for  small  arguments  is 

Y,(x)  - - (B-29) 

wx" 

giving 

V ■ -jTf,  tY,.,  --2X1  (B-30) 

v'  ♦ ■ J)  ■ '‘i  * 1 - r ' • irrT  (B-31) 


Substituting  these  relations  into  Eq.  (B-12) 

where  Yu  > Yn(ap,i).  Subatitution  of  this  expression  in  Eq.  (B-13)  gives  i in  the  case 
^ B 0.  the  term  aPn  becomes  afl/c,  and  one  obtains 


afc,  * e*^  ^ ~ 


Eq.  (B-16)  gives 
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The  vital  condition  (B-26)  for  the  validity  of  above  results  becomes 

1*6,1  (B-35) 

If  “ 0 or  (n  ~ 0,  then  p,,  » n/c  and 

|a6i|«(n  - 1)  ^ (4s  -0)  (B-35a) 

This  condition  now  shows  clearly  that  the  results  of  this  appendix  are  valid  for  small  4| 
where  they  are  to  be  applied.  The  second  condition  (B-28)  does  not  simplify  for  small 
arguments. 

Conclusion; 

Substituting  Eq.  (B-27)  into  Eq.  (32)  and  noting  that  A is  not  a function  of  z,  the 
displacement  can  be  expressed  in  terms  of  the  displacement  w(0)  at  z « 0, 

it|t 

w(z)  = w(0)e  ' (B-36) 

where  the  complex  q\iantity||  follows  from  Eqs.  (B-12,  13).  The  relation  (B-36)  is  valid 
only  if  the  condition  (B-26)  is  strictly  satisfied,  and  Eq.  (B-28)  is  at  least  not  badly  vio- 
lated. For  small  arguments  Eqs.  (B-32,  33  and  35)  can  be  used  instead  of  Eqs.  (B-12,  13 
and  26),  respectively. 
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